Rapid State-Reduction of Quantum Systems Using Feedback Control 
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We consider using Hamiltonian feedback control to increase the speed at which a continuous 
measurement purifies (reduces) the state of a quantum system, and thus to increase the speed of 
the preparation of pure states. For a measurement of an observable with A'' equispaced eigenvalues, 
we show that there exists a feedback algorithm which will speed up the rate of state-reduction by 
at least a factor of 2{N + l)/3. 
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To perform information processing using quantum sys- 
tems the systems in question must be prepared in pure 
states . Due to environmental noise such systems often 
exist naturally in mixed states, and either a process of 
cooling or measurement must be used to purify them. 
While it is often assumed for the sake of simplicity that 
measurements arc instantaneous, this is not the case in 
reality. All physical processes happen on some finite time 
scale, and measurement is no exception. Thus every von 
Neumann measurement is in fact a continuous measure- 
ment which has been run for sufficiently long. The speed 
at which one can measure a system will therefore place 
limits on the speed at which measurement can be used 
as a tool to prepare initial states. 

Given a fixed measurement rate, it is possible to use 
Hamiltonian feedback [sl Ij, lu Li la| during the mea- 
surement to increase the speed at which the system is 
purified. For a single qubit this problem was analyzed 
by Jacobs 0. For this case the optimal feedback algo- 
rithm will produce a factor of 2 increase in the speed 
at which a given level of purity can be obtained. This 
speed-up factor is achieved in the limit in which the de- 
sired purity is high, being the limit in which the measure- 
ment time is large compared to the measurement rate. 
Shorter measurement times result in a smaller speed-up. 
The speed-up is a quantum mechanical effect, in that the 
equivalent measurement on a classical bit cannot be en- 
hanced in this manner; for a two-state system a speed-up 
is only possible if the system can exist in a superposition 
of the eigenstates of the measured observable during the 
process. Specifically, the optimal algorithm consists of 
measuring a^, and using Hamiltonian feedback to con- 
tinually rotate the state of the system as needed so that 
it remains diagonal in the <Tx basis throughout the mea- 
surement. This requires a process of real-time feedback, 
since the required stream of unitary operations depends 
upon the continuous stream of measurement results. 

Here we consider using feedback to enhance a mea- 
surement on a system of arbitrary size. This problem is 
considerably more complex than that of a single qubit, 
because one cannot directly derive the analytical results 
which are available in that case We will consider 



a continuous measurement of an observable with N dis- 
tinct, equally spaced eigenvalues (i.e. one that distin- 
guishes the A^-states), present a feedback algorithm, and 
prove a lower bound on its performance. Using this 
lower bound we show that the speed-up factor which 
the algorithm provides in the long time limit is at least 
2(iV-|-l)/3. Thus for large systems Hamiltonian feedback 
can greatly enhance the rate of state-reduction. 

The evolution of a system under the continuous mea- 
surement of an observable X is given by 0, 0, 0, Q, 



dp = -7[X, [X,p]]dt+^2j{Xp + pX - 2{X)p)dW, (1) 

where p is the density matrix of the system and dW is 
an increment of the Wiener process, describing driving 
by Gaussian white noise, and 7 is a positive constant 
often referred to as the measurement strength. The con- 
tinuous stream of measurement results, y(t), is given by 
dy(t) = {X(t))dt + dW/^/Wj. The effect of such a mea- 
surement in the long time limit is to project the system 
onto one of the eigenstates of X, producing a von Neu- 
mann measurement in the basis of X. 

The rate at which the measurement distinguishes be- 
tween two eigenstates of X is proportional to 7 multiplied 
by the difference between the respective eigenvalues. In 
order to able to purify any initial state the observable 
must therefore possess non-degenerate eigenvalues. As 
mentioned above we will also assume that the eigenval- 
ues of X are equally spaced. This is true, for example, if 
the system has total angular momentum j, and the ob- 
servable is the z-component of angular momentum, Jz. 
In this case there are A^ = 2j + 1 states, with eigenval- 
ues Further, since Eq.J^I is invariant under 
the transformation X X' = X -\- al for any real a, 
an overall shift in the eigenvalues of the measured ob- 
servable has no effect on the measurement. As a result, 
without loss of generality we may always assume that X 
is traceless. Under the restriction that the eigenvalues 
are equally spaced, we may therefore always assume that 
X is proportional to Jz, where j = {N — l)/2. While 
the constant of proportionality is unimportant for the 
purposes of our analysis, we will define X so that its 
eigenvalues are separated by 1/A^. Thus X = Jz/N. 
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Using the algorithm in |9| as a guide, we will consider 
a feedback algorithm consisting of a measurement of X, 
in which the feedback is used to rotate the state of the 
system so that the eigenbasis of the density matrix re- 
mains unbiased with respect to the eigenbasis of X dur- 
ing the measurement. Recall that two bases are unbiased 
with respect to one another when the moduli of the ele- 
ments of the matrix which transforms between them are 
all identical . We will refer to the eigenbasis of X as 
the measurement basis. To implement such an algorithm 
one applies a Hamiltonian in each time interval to rotate 
the state of the system so as to cancel the infinitesimal 
change in the eigenbasis caused by the measurement in 
that time interval. In analyzing the performance of this 
algorithm, we will use the impurity, L = 1 — Tr[p^] [2(1 |. 
as a our measure of mixedness, rather than the von Neu- 
mann entropy, because it has a very simple analytical 
form. In addition, the impurity is equal to the von Neu- 
mann entropy in the limit of high purity. 

To begin with we calculate the instantaneous derivative 
of L when the measurement basis is unbiased with respect 
to p. In this case, if we work in the eigenbasis of p, all 
the diagonal elements of X are the same, and equal to 
(X), and we have further that Tr[X"p™] = (X")Tr[p™]. 
Using this fact in Eq.QJ we find that the evolution of L 
is particularly simple: 

dL = -8-fTr[XpXp]dt. (2) 

That is, when the measurement basis is unbiased with 
respect to p, the evolution of the impurity over the suc- 
ceeding infinitesimal time step is deterministic. In the 
absence of any feedback, the action of the measurement 
changes the eigenbasis of p, so that the measurement docs 
not remain unbiased. Under our feedback algorithm, 
however, Eq.Q remains true, and the purity therefore 
evolves deterministically. 

The key to obtaining a lower bound on the performance 
of our feedback algorithm is the following lemma. 

Lemma 1. Consider a traceless observable X of an N 
dimensional quantum system, whose eigenbasis is unbi- 
ased with respect to a density operator p, and consider 
the N\ unitary operators whose action is to permute the 
eigenvectors of p. There exists a set of N\ operators 
X^™) = UmXU^ all unbiased with respect to the state p 
of a system, such that when these operators are measured 
simultaneously with strength 7, the instantaneous deriva- 
tive of L = 1 — Tr[p^] is given by dL/dt = —SkjNlL, 
where k = Tr[X^]/{N{N - 1)). 

Proof. Consider a simultaneous measurement of M trace- 
less operators X^™^ unbiased with respect to p. From 



Eq.(|21l, the derivative of the impurity is 

dL = -87d<^rr[X(™V^^"V] (3) 

rn 

m ij 

where the xj"^^ are the matrix elements of the X^"^'> 
''J 

in the eigenbasis of p, and the Ai are the eigenvalues 
of p. Since the X'™) are traceless, xj^l"'' = 0, Vi, m. 
Now consider that the X^™^ are the iV! permutations of 
a single unbiased traceless operator X. In this case, for 
^^Jwe have l^if^l' = ^ 2)! Y.i^, = 
{N — 2)! Ej,; l^fcjp = c. Naturally when i = j we have 
E™l^ir^P = 0. As a result 

dL = --8-/dt^c{l~ Sij)X,Xj (5) 

= -8-fcdtL. (6) 

Permuting the elements of the operator X when writ- 
ten in the eigenbasis of p corresponds to a basis change 
which involves permuting the vectors of the eigenbasis of 
p. Such a transformation leaves X unbiased with respect 
to p, and thus all the iV! operators X*^™' arc unbiased 
with respect to p, and related to X by unitary transfor- 
mations. 

One can also calculate the constant c by using 
X;„Tr[X(™)pX('")p] = cdtL and setting p = I/N. The 
result is c = (iV - 2)!Tr[X2]. □ 

This result shows us two things. The first is that if 
we measure simultaneously A^! unbiased rotated versions 
of X , and combine this with a feedback algorithm which 
maintains the eigenvectors of p, the impurity will decay 
as 

L{t) = e-8^(^-2)'T'-[^Y=]t^(Q)^ 

The second is that, because we know that the derivative 
of the impurity when summed over all of the N\ mea- 
surements is —SjNlkL, we know that there is at least 
one term in this sum for which the derivative is greater 
than or equal to —8^kL. That is, given an unbiased op- 
erator X, and a density matrix p, we can always find a 
transformation T (which permutes the eigenbasis of p), 
such that measuring TXT'^ will provide an instantaneous 
change in the impurity greater than or equal to —8jkL. 
Of course, we can equivalently apply the transformation 
to p, rather than to X, which simply involves permuting 
the eigenvalues of p. 

We can now return to our feedback algorithm. Recall 
that this involves measuring a single unbiased operator 
A, and applying a unitary operation to the system at 
each time step to keep the eigenvectors of p unchanged. 
However, using our unitary operation we can also per- 
mute the eigenvectors of p. Therefore, at each time step, 
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we can choose the permutation T so as to maximize 
the derivative of the impurity. As we have now shown 
that this maximum derivative is bounded from below by 
—8jkL, we see that our feedback algorithm will produce 
an impurity Lfb(t) such that 



Libit) < e 



-87Ti-[X- 



i*/(^'-^)l(o). 



(8) 



We now have a lower bound on the performance of 
our feedback algorithm. However, what we really wish 
to know is the enhancement that the algorithm provides 
over the unaided measurement. In this case the evolu- 
tion is more complex than that under the feedback al- 
gorithm because it is stochastic. However, Eq.([3l can 
be solved using, for example, the technique developed in 
reference [i^ ■ If the state of the system is initially com- 
pletely mixed, then the solution can be written, up to a 
normalization constant, as 



p{t, v) oc e' 



-4'rt{X^+2vX) 



(9) 



Here v is the measurement record integrated up until 
time t, and divided by t. That is. 



v{t)^^ f dy = \ j\x{s))ds + J 



r* dw{s) 



(10) 



As will become clear shortly, v{t) is the final result of the 
measurement after measuring for a time t: if v is close n 
at time t, and t 3> A^^/(87), then the system has been 
very nearly projected onto the eigenstate of X with eigen- 
value n (recall that n takes the values —j/N, . . . ,j/N). 
The probability density for v is 
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(11) 



Thus for times long compared to I/7, v is sharply peaked 
about the N values —j/N,...,j/N, each peak having 
width l/^/Wyi. 

The density matrix remains diagonal in the eigenbasis 
of X, and its diagonal elements are 



Pn{t) 



-4-yt(v—nY 



(12) 



where the normalization \s M = e"^'''*^^"") . We 
see that the size of the n*^ element is determined by the 
proximity of n to u. When v is close to n, and t is large, 
Pn is much greater than all the other elements. 

We wish to calculate the average impurity in the 
long time limit, so as to compare it with the impu- 
rity achieved under the feedback algorithm. That is, 
(P(t)) = 1 - (Tr[p2]) = 1 _ /^Tr[p2]p(„)d„. We will 

show first that for long times the average purity, (Tr[p^]) 
for an N state system is bounded from above by that for 
a measurement of az/'i-N on a two-state system. Note 
first that for long times P{v) is sharply peaked at the 



N eigenvalues of X, so that (Tr[p^]) can be written as 
the average of N integrals, in each of which we only need 
consider values of v which are close to a single value of 
n. Since this is true for every N , the only difference be- 
tween the average purity for different values of N comes 
from the purity Tr[p^]. Now consider the two integrals 
corresponding to the extreme values oi n {n = j/N and 
n = —j/N). We note that for any value of N, if we 
truncate the density matrix so that we keep only the 
two largest elements (and renormalize) , then the value of 
Tr[/9^] in these integrals becomes precisely that describ- 
ing a measurement of az/'2N on a two state system. The 
factor of 1/iV stems from the fact that adjacent eigen- 
values of X are separated by 1/N. Further, the process 
of truncation increases the purity. This follows imme- 
diately from the theory of majorization applied to the 
diagonal elements of p, and the fact that the purity is a 
symmetric convex function of these elements [23 • Thus, 
the two integrals corresponding to the extreme values of 
n are bounded from above by those for a measurement of 
az/2N. Further, it is immediately clear that for every N, 
the value of Tr[/9^] under the (A^— 2) integrals correspond- 
ing to the non-extreme values of n, upon truncation, give 
a smaller purity than those for the extreme values (this 
is because for these values of n, there are adjacent values 
on both sides, and not merely one side). Since the total 
average purity is simply the average value of the A^ inte- 
grals, for every A^ this average is bounded from above by 
that for a measurement of az/2N. Thus for long times 
the value of L for such a measurement provides a lower 
bound on a measurement of X for N > 2. (In fact, nu- 
merical results show that this is true at all times - as 
indeed one would intuitively expect - and not merely in 
the limit in which t > A'^/(87).) 

The average value of the impurity for a measurement 
of az/2N may be written as 
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+ 00 



-xV(2t) 



cosh{y^x/N) 



dx. 



(13) 



Noting that the integral in this expression is independent 
of t for t > Af2/(87), we have, in this hmit, P2(i) = 
e-T*/^'c/(^/8^), where C = 1 / cosh{y^x/ N)dx. 

Now that we have an upper bound on L under our 
feedback algorithm, and a lower bound on this for the 
unassisted measurement, we can derive a lower bound 
on the speed-up factor, 5*, provided by the algorithm. 
Specifically, 5* is the ratio between the time taken to 
achieve a given target value of the purity under the un- 
aided measurement, <ni, and that required when using 
the algorithm, tfi,. The result is 



1 1 

S - 8iV2fc 



(14) 



where C = 
Tt[X^]/{N{N 



V8^{N - 1)/(CA^). 
1)), evaluating Tr[X^] 
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Tr[j2]/iv2 
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Target Impurity 



FIG. 1: We plot here a lower bound on the speed-up factor 
provided by our feedback algorithm as a function of the target 
impurity, when the initial state is completely mixed. The solid 
line gives the result for system dimension N = 2, the dashed 
line for = 3 and the dot-dashed line for = 4. 

Y.i=-jn^ = iN + l){N--l)/{12N), and taking the limit 
of large t, we obtain 

5>^(iV + l). (15) 

This expression is valid in the limit of high target 
purity. We have also calculated numerically the lower 
bound on the speedup factor as a function of the target 
purity to examine how it approaches the above limit, and 
we plot the results for iV = 2, 3 and 4 in figure We find 
that this lower bound increases monotonically towards its 
limiting value as the target purity is increased. 

Before we finish we describe in a little more detail the 
feedback Hamiltonian involved in implementing the al- 
gorithm. To implement the feedback algorithm the con- 
troller is required to apply a unitary operation to the 
system so as to preserve the basis in which the density 
matrix, p, is diagonal, and, if necessary, apply a uni- 
tary operation to reorder the eigenvalues of p. To cal- 
culate the first unitary the controller diagonalizcs dp. 
Since dp has two terms, one proportional to dt and 
the other to dW, this unitary will in general take the 
form U{t, dt) = eKHi{t}dt+H2{t)dW) ^^^^ operators Hi 
and i?2- The feedback Hamiltonian is thus, in general, 
H{t) = Hi{t) + H2{t){dW/dt), and therefore contains a 
term proportional to the noise part of the measurement 
record y{t). This kind of quantum feedback has been 
studied extensively by Wiseman and Milburn The 
second unitary, that required to reorder the eigenvalues 
of p should ideally perform the operation in a single time 
step. This requires a Hamiltonian proportional to 1/dt, 
and since dt should ideally be small on the time scale of 
the dynamics due to the measurement, this would require 
a Hamiltonian large compared to 7. 



A number of interesting open questions remain. While 
we have found a lower bound on the performance of our 
algorithm, its actual performance may be higher. In ad- 
dition, we do not know whether this algorithm is op- 
timal. Further, in the above analysis we have focused 
on the ability of feedback to increase the speed of state- 
reduction under the assumption that there is no signifi- 
cant limitation on the feedback Hamiltonian. It will be 
interesting to consider the performance of the algorithm 
when specific limits are placed on the feedback Hamilto- 
nian, and this will be the subject of future work. 
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